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There is evidence for an excess of gamma rays with O(GeV) energy coming from the Galactic 
Center in data from the Fermi Telescope. The spectrum of the excess is well fit by 30 GeV dark 
matter annihilating into a pair of b quarks with a cross section expected from a thermal relic. For an 
explanation of this excess, we study a renormalizable model where dark matter couples to the SM 
via a pseudoscalar that mixes with the CP-odd Higgs boson of a Two Higgs Doublet Model. We 
report the constraints on this model from direct detection, Higgs decays, and rare B meson decays. 


I. INTRODUCTION 

Eventhough dark matter (DM) accounts for more than a quarter of the energy density of the Universe, the 
only evidence of its existence comes from gravitational observations such as galactic rotation curves, cluster 
mergers and the cosmic microwave background (see, e.g. [1] and references therein). There are three ways 
to search for non-gravitational interactions of DM: (i) Collider searches where DM is produced through SM 
interactions, (ii) direct detection, and (iii) indirect detection where products of DM annihilation or decay are 
searched for. 

There is evidence for indirect detection of DM annihilaition, several groups finding an excess of gamma rays of 
energy ~l-3 GeV in the region of the Galactic Center observed by the Fermi Gamma Ray Space Telescope mm- 
Astrophysical sources are unlikely explanations for this excess due to its spectrum and morphology [1]. On the 
other hand the spectrum of the excess has been fit by DM annihilating to a number of final states, depending 
on its mass, notably 10 GeV DM annihilating to t+t“ (and possibly other leptons) and 30 GeV DM to bb 
[5 . The size of the excess is compatible with an annihilation cross section roughly equal to that expected for a 
thermal relic, {avre\) = 3 x 10“^® cm^/s, suggesting that it is actually the result of DM annihilation. Here we 
focus on a 30 GeV DM that annihlates to b quarks. The effective pseudoscalar operator 

( 1 ) 

where y is the DM, has been the most studied interaction to explain the excess [6] since its direct detection 
cross section is velocity suppressed and spin dependent, rendering it safe from direct detection experiments. 

In the rest of this paper we summarize the results of [7], as presented in Toyama International Workshop 
on Higgs as a Probe of New Physics 2015, constructing a UV complete model that produces the pseudoscalar 
interaction in Eq. [l] where DM-SM interactions are mediated by a pseudoscalar that mixes with the CP-odd 
Higgs boson of a Two Higgs Doublet Model (2HDM). We also report the constraints on this model. 


II. THE MODEL 

To construct a pseudoscalar-pseudoscalar interaction between DM and the SM, we expand the Higgs sector of 
the SM with a second Higgs doublet, which has enough degrees of freedom to mix with a pseudoscalar mediator. 
We take the DM to be a Dirac fermion, y, with mass coupled to a real, gauge singlet, pseudoscalar mediator, 
Gq, through 


Cdark = VxO-oXh^X- 


( 2 ) 


The mediator couples to the SM via the Higgs portal in the scalar potential which is 

V = k 2 HDM + 2 ^ao®o A iBaoH\H 2 + h.c. (3) 

with Hi ^2 the two Higgs doublets. We assume that £dark and V are GP-conserving (i.e. B and y^ are both real, 
and there is no CP violation in V 2 hdm)- In this case, qq does not develop a VEV. We write the most general 
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CP-conserving 2HDM potential as 

.,2^ 2 


V2HDM — Ai ( H\Hi -] + A2 ( H2H2 -) + A: 


hIh, - ^ ) + ( H^H2 - ^ 


rt i 


( 4 ) 


+ A 4 


h\h^ 


(hIh^) - (h\H2) 


+ As 


Re [hIh.] - 


V 1 V 2 

2 J 


+ Afi 


Im 


(h\H2) 


with all Ai real. We have also imposed a Z 2 symmetry under which Hi —>■ Hi and H 2 —> —H 2 to suppress 
flavor-changing neutral currents, which is only softly broken by V 2 HDM and Vport- Tke potential is minimized 
at {Hi) = Vil\f2, i = 1, 2, and the W and Z masses fix v\Pv\ = v^ = (246 GeV)^. The angle P is defined by 
tan/3 = V 2 1Vi- In unitary gauge we can decompose the doublets as 


•■=-:=( V 

y/2 \ Vz + p^ + ir]i J 


( 5 ) 


The spectrum contains a charged Higgs (H^), two GP-even Higgses {H,h), and two GP-odd states {A, a) that 
are a mixture of the GP-odd Higgs Aq and the pseudoscalar mediator oq, 

H^ = sin j3 (j)f— cos fS (f>^, H = cos a pi -I- sin a p 2 , h = — sin a pi -I- cos a p 2 , 

H = cos6*Ho + sin6*ao, a = — sin0Ho + cos0 oq, (6) 


where Aq = sin /3 — cos /3 772 - For details of the scalar sector, including definitions of the mixing angles, see [7]. 

The DM coupling to the mediator in Eq. ([^ is expressed in terms of CP-odd mass eigenstates. 


^dark = Vx (cos Oa + sinOA) x*7 X- 


( 7 ) 


We will work in a Type H 2HDM, where H 2 couples to up-type quarks and Hi couples down-type quarks and 
leptons. The couplings of the neutral scalar mass eigenstates are then rescaled from the SM Higgs values by 


th 

Se,d 


bill UL ^ CUb Ct ^ 

COSO ^ sin a ^ 


tan P COS 6, ^ = — tan /3 sin 9, 

cot /3 cos 9, = — cot P sin 9. 


( 8 ) 

(9) 


To simplify the analysis, we work close to the decoupling and alignment limit of the 2HDM where a ~ P — 'kI2 
and m/i <C mn — rn[j± ~ mAg- Since h has SM-like couplings in this limit, we identify it with the 125 GeV 
Higgs boson. 


III. DARK MATTER ANNIHILATION 


For rria <C nriA, the dark matter annihilates to SM particles primarily through s-channel a exchange. The 
velocity averaged annihilation cross section for xx 8M in the nonrelativistic limit is 


y^sIgUn^pm^ 

CTUrel) = o 4 1 

Stt 

( 10 ) 


4to^ 


w 

a 

2 


m; 


f=b,T,... 


Nc- 


m% 


'1 -( ~ 3 X 10"^® 


cm 


y^z sin 29 tan /3 ^" 
2.4 J 


for the experimentally favored DM mass of 30 GeV, taking ma = 100 GeV and ignoring Fa. 

It is possible to achieve values of the annihilation cross section compatible with the gamma ray excess and 
relic density constraints with modest values of the mixing angle 9, provided tan/3 is somewhat large. At this 
value of my_ and for tan P larger than a few, the bb final state accounts for about 90% of the annihilation cross 
section with making up nearly all the rest. This is in line with what is suggested by fits to the gamma 

ray excess. 
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FIG. 1: Regions of mixing angle 6 vs. ma that are ruled ont or suggested by various measurements. We have fixed 
rrifj ~ rriA = 800 GeV, tan/3 = 40, a = /3 — tt/2, and = 0.5. The area that gives an annihilation cross section of 
{aVrei) = 1 — 5 X 10“^® cm®/s as indicated by fits to the gamma ray excess is between the solid black lines (shaded in 
green). Other shaded regions are (projected) exclnsion regions coming from various current (future) measurements. See 
text and also [7] for details. 


IV. CONSTRAINTS ON THE DARK SECTOR 


In this section we investigate the limits on the mediator mass and the mixing angle between the mediator and 
the pseudoscalar of the 2HDM. We use the benchmark parameteres mn = mH± — tua = 800 GeV, tan/3 = 40, 
a = /3 — 7r/2, and = 0.5. We consider constraints from direct detection, Higgs decays and rare H-meson 
decays. As discussed in [7], this model is not well constrained by monojet searches due to the suppressed 
coupling to top at large tan /3. Our main results are summarized in Fig. [l] 

1. Direct Detection 

Pseudoscalar exchange between y and quarks leads to (highly suppressed) spin-dependent scattering of 
the DM on nuclei. However, spin-independent interactions are generated at one-loop. For our benchmark 
values the dominant contribution to this comes from Higgs exchange triangle diagrams, leading to a 
spin-independent direct detection cross section of 


crsi ~ 2.2 X 10 "^® 


2 ( V / 50 GeV y / TOx 

VSOOGeVy \ ma J V30 GeV/ ^0.1/ V0.5/ 


( 11 ) 


In Fig. we show the area of parameter space ruled out by the LUX limit of 8 x 10 cm^ [H]. We also 
show the area that can be probed by a future cross section limit of 10“^® cm^. 


2. Higgs Decays 

If ma < mh/2, Higgs can decay into a pair of a’s with the decay rate 


F (h —)■ aa) 


{m\ — ma) ^ sin'* 20 
32TTmhv‘^ 



~ 840 MeV 


mA 


.800 GeV 



( 12 ) 


for ma mh,A- This can impact LHG measurements of h, with a SM width of 4 MeV, which is SM-like 
to 10 - 20%. ’ 

For ma < TOh/2 ~ 2m^, the pseudoscalars will primarily go to b quarks with a small branching to r and 
p pairs. The h ^ aa ^ Ab,2b2p signal will contribute to h ^ bb searches [9]. A GMS search in WjZ- 
associated production at 7 and 8 TeV, pp —>■ {WjZ) + {h ^ bb) [TU], sets a limit Hr (h ^ aa ^ 46) < 0.7 
for 2mb < ma < 15 GeV, and Hr {h ^ aa ^ 2b2p) < 10“^ for ma > 25 GeV. 
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Since we are in the decoupling limit, the production cross section of the Higgs is unchanged from its SM 
value in this model. Therefore there are strong limits on unobserved final states, such as aa, that would 
dilute the signal strength in the observed channels. Given current data, this limits Br {h —>■ aa) < 0.22 [llj . 

We show in Fig. the limits coming from Br (h —>■ aa —>■ 46) < 0.7 as well as the indirect constraint 
Br {h —>■ aa) < 0.22 and the limit that can be set by a future measurement of Br (/i —?> aa —>■ 262/i) < 10“"^. 

3. B Physics Constraints 

A light a can also potentially be constrained by its contributions to the decay Bg —)■ ■ For TOq «C mz, 

the correction due to s-channel a exchange can be simply written as [12] 


Br {Bs ~ Br {Bg )gj^ X 


mbmB„tlsl f {xt,yt,r) 
- "la F (xt) 


(13) 


with Xt = mf/rn^, yt = r = and / and Y are loop functions that are defined in 

[7]. The average of the LHCb and CMS measurements of this mode is Br {Bg —>■ yYpT) = (2.9 ± 0.7) x 
10“® [13]. This should be compared against the SM prediction, which we take to be (3.65 ± 0.23) x 
10“® [13]. This offers a strong test of the model, especially for a light a, which we show in Fig.[^ 


V. CONCLUSIONS 

In order to explain the gamma ray excess from the Galactic Center, we have studied a 2HDM where the 
pseudoscalar mediator mixes with the CP-odd Higgs, giving rise to interactions between DM and the SM. 

A spin-independent cross section for direct detection arises at one-loop level in this model and is well below 
the current bound of 8 x 10“^®cm^ at a dark matter mass of 30 GeV. We also consider decays of the 125 GeV 
SM-like Higgs boson involving the mediator. If the mediator is light h —^ aa —>■ 46, 262/i can be constraining 
with data from the 14 TeV LHC. Additional contributions to Bg —>■ y~ in this model eliminate some of the 
favored parameter space for ma <10 GeV. 
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